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We study the crossover from Ising-like to classical critical behavior as a function of the range 
R of interactions. The power-law dependence on R of several critical amplitudes is calculated from 
renormalization theory. The results confirm the predictions of Mon and Binder, which were obtained 
from phenomenological scaling arguments. In addition, we calculate the range dependence of several 
corrections to scaling. We have tested the results in Monte Carlo simulations of two-dimensional 
systems with an extended range of interaction. An efficient Monte Carlo algorithm enabled us to 
carry out simulations for sufficiently large values of R, so that the theoretical predictions could 
actually be observed. 
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I. INTRODUCTION 



> 

Q\ ■ As is well known, the critical behavior of a physical system strongly depends on the range of the interactions. The 
longer the range, the stronger critical fluctuations will be suppressed. In the limit of infinite range we recover classical 
or mean-field-like critical behavior. For interactions with a finite range, however, fluctuations remain very important 
and essentially modify the critical behavior. As follows from the Ginzburg criterion jlj, sufficiently close to the critical 
temperature T c nonclassical critical exponents apply for any finite interaction range R. This crucial difference between 
finite and infinite R implies a crossover from one type of critical behavior to another as a function of R. Such crossover 
phenomena are of great interest for a wealth of critical phenomena. They occur, e.g., in polymer mixtures (see Ref. 
and references therein), as a function of the chain length, and gas-liquid transitions, as a function of the difference 
between the temperature and the critical temperature. The explanation of these phenomena in terms of competing 
fixed points of a renormalization transformation is one of the important features of the renormalization theory (see, 
e.g., Ref. |). 

In Ref. [ 1 1 , Mon and Binder have already studied crossover as a function of R within the context of finite-size scaling, 
motivated by the crossover in polymer mixtures. They predicted that the critical amplitudes of scaling functions 
display a singular dependence on R. The various power-law dependencies were obtained from phenomenological 
crossover scaling arguments. In this paper, we will derive this dependence on R from a renormalization description 
of the crossover from classical to nonclassical critical behavior. The first part of the renormalization trajectory is 
governed by the Gaussian fixed point, which is unstable for d < 4. The corresponding scaling relations have been 
derived by Rikvold et al. S. Sufficiently close to criticality, the final part of the renormalization trajectory is governed 
by the Ising fixed point. The resulting relations are in complete agreement with the predictions from Ref. Q. In 
addition, we obtain the R dependence of the leading corrections to scaling and derive from renormalization arguments 
a logarithmic factor in the shift of the critical temperature. This factor was already conjectured in Ref. (J]. 

It is interesting to note that the physical mechanism leading to the singular range dependence of scaling functions 
is closely related to that leading to the violation of hyperscaling for d > 4. The latter effect is caused by a singular 
dependence of thermodynamic quantities on the coefficient u of the cf> 4 term for u — ► in a Landau-Ginzburg- Wilson 
Hamiltonian. In other words, u is a so-called dangerous irrelevant variable (see, e.g., Ref. [^) for a more detailed 
discussion). Here, as we will see, the fact that this coefficient becomes small for large values of R plays again an 
essential role, although u is relevant for d < 4. 

Furthermore, we present new Monte Carlo results for two-dimensional Ising models with an extended range of 
interaction. A serious problem associated with such simulations is that the simulation time tends to increase rapidly 
with the number of interacting spins. However, a large interaction range is crucial to observe the predicted R 
dependencies, as will follow from the renormalization description. The maximum range that could be accessed in 
Ref. jl| was too small to verify the theoretical predictions. We overcome this limitation by means of an efficient 
cluster Monte Carlo algorithm 0, in which the simulation time per spin is practically independent of the range of 
the interactions. 
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The outline of this paper is as follows. In section Q we derive the R dependence of critical am plitu des from 
renormalization theory. These results are verified by Monte Carlo simulations, presented in section III. We end 
with our conclusion in section IV. Two technical issues, namely the Fourier transform of a spherical distribution 
of interactions in a general number of dimensions and the application of the cluster algorithm to medium-range 
interactions, are addressed in Appendices [A| and [b|, respectively. 



II. RENORMALIZATION DERIVATION OF THE DEPENDENCE OF CRITICAL AMPLITUDES ON 

THE INTERACTION RANGE 

A model with long-range interactions which has attracted much attention is that in which the spin-spin interactions 
decay algebraically as a function of the distance r between the spins, K(r) = Ar- d - a (a > 0), where d is the 
dimensionality and A > the interaction strength (see, e.g., Refs. |g|,||j6|). For d < 4, this model displays an 
interesting continuous variation of critical behavior as function of a: for < a < d/2 the critical behavior is classical 
(mean- field-like), for d/2 < a < 2 — rj sr the critical exponents vary continuously, and for 2 — r/ sl - < a, the interactions 
decay fast enough to yield short-range (Ising-like) critical behavior. Here, ?7 sr denotes the exponent rj in the d- 
dimensional model with short-range interactions. The algebraical decay of the interactions is responsible for the 
existence of an intermediate regime between Ising-like and classical critical behavior. In this paper we focus on a 
different way to interpolate between the long-range (mean-field) limit and short-range models. Instead, we choose 
ferromagnetic interactions which are constant within a range R and zero beyond this range. Thus, we have the 
following Hamiltonian 



H/k B T= -J2J2 Kd ^ 



(i) 



where the spin-spin interaction -Kd(r) = cR~ d for |r| < R and the sums run over all spins in the system. This 
Hamiltonian displays physical behavior that is different from the power-law case. In particular, the intermediate 
regime with variable exponents is absent, and mean-field critical behavior is restricted to the infinite-range limit. We 
analyze the influence of the range R within the context of renormalization theory, starting from a generalized Landau- 
Ginzburg- Wilson (LGW) Hamiltonian, where the [V(/>(r)] 2 term normally representing the (short-range) interactions 
is replaced by an interaction term with spin-spin coupling ( |Al| ) , 



dr ■ 



dr' 



-r'\<R 



(r) + i V 2 (r) +u o 4 (r) 



(2) 



As a consequence of the normalization factor R~ d , the critical value of the temperature parameter c depends only 
weakly on R. The first integral runs over the volume V which contains TV particles. We adopt periodic boundary 
conditions. The Fourier transform of the interaction is calculated in Appendix |Aj. It leads to the following momentum- 
space representation of the Hamiltonian 



2n 



d/2 



— YYY> 

4N t—' ^ ^ 



'k, '.'k»"k3«>-ki-k2-k3 — ho\l "^k=o 



(3) 



J„ is a Bessel function of the first kind of order v. The wavevectors are discrete because of the periodic boundary 
conditions. Furthermore, we restrict the wavevectors to lie within the first Brillouin zone, which is reminiscent of the 
underlying lattice structure. The interaction term can be expanded in a Taylor series containing only even terms in 
kR. This means that we will be mainly concerned with the term of order (fci?) 2 , because higher-order terms will turn 
out to be irrelevant. The constant term in the Taylor series is absorbed in v and the coefficient of the quadratic term 
as a factor in c. This yields a new Hamiltonian 



W t (0k)/fcBT = ^ Y, [ cR2,{2 + v ] <^k0-k + ^ Y Y Y <ftk x <^k 2 4>k 3 <ft-ki-k 2 -k 3 - h J~^4>k= o 



(4) 



ki k 2 k 3 



Since we are free to choose the scale on which the fluctuations of the order parameter are measured, we may rescale 
(p — ¥ ip = \ftR(\>. This is generally the most convenient choice because the dominant fc-dependent term becomes 



2 



independent of This leads to 

n(^)/k B T = ~ x [ fc2 + ^2] V'k^-k + 4 _ 2 ^ 4jV XI X 51 V'k 1 V'k 2 V'k3V'-k 1 -k 2 -k3 - -fey^v-k^o . (5) 

k ki k 2 k 3 V 

The parameter c is merely a constant, independent of the range, and in order not to be hampered by it in the future 
analysis, we absorb the various powers of it in ro = v/c, u = uq/c 2 , and h= ho/ \/S. Now, ro assumes the role of the 
temperature parameter. 

If the range i? is large, the coefficient of the tp 4 term is relatively small and hence the critical behavior of the system 
is determined by the Gaussian fixed point. Under a renormalization transformation with a rescaling parameter I the 
Hamiltonian thus transforms as 

H'(^,)/k B T> = lj2[ k ' 2 + "fA Hr#-u + ^^EEE<l^^^- k i- k ; 

k' k'j k 2 k 3 

-^l 1+d/2 ^o ■ (6) 

Here i/j' k , = l -1 ^, k' = ki, the sums run again over the full Brillouin zone, and N' = Nl~ d . For d < 4, the ip 4 term 
grows and the system moves away from the Gaussian fixed point /Zq (see Fig. [l]). At present, we are interested in 
the flow from the neighborhood of the Gaussian fixed point to that of the Ising fixed point. Thus we remain close 
to the critical line connecting the two fixed points and the temperature field parametrized by ro remains small. The 
crossover to Ising-like critical behavior occurs when the coefficient of the ip 4 term is of the same order as that of the 
k 2 ip 2 term, which is unity, i.e., when I = Iq = R 4 ^ 4 ~ d \ We shall refer to Iq as the crossover scale. 

By comparing the coefficient of the ip 4 term to that of the roip 2 term, it is possible to derive a criterion that states 
for which temperatures the critical behavior will be Ising-like and for which temperatures it will be classical. This 
is the well-known Ginzburg criterion El, which can also be derived from Eq. (||) (see, e.g., Ref. @ p. 107]). One 
expects the Gaussian fixed point to dominate the renormalization flow if, irrespective of I, the ip 4 coefficient is small 
compared to the temperature coefficient. Thus, one requires the scaled combination uR~ 4 l 4 ~ d / (roR~ 2 l 2 )( 4 ~ d ^ 2 to 

be small, or r (4 ~ d)/2 R 1 *^ 1 > 1 (cf. also Ref. g Eq. (3)]). 

Since we are now in the neighborhood of the Ising fixed point, we continue renormalizing our Hamiltonian with 
nonclassical renormalization exponents yt, y,h, and y x . To leading order, it will transform as follows, where b denotes 
the rescaling factor of our new transformation, 



n"(K„)/k B T" = \ £ [k" 2 + R 2d ^- d \ b ytf + rg)] + b -^- EEE ^k^-K-K-K 

k" k '/ k' 2 ' k' 3 ' 

-hR* d /^\J^by^,, =Q . (?) 

We have introduced the coefficients fo and u, which denote the location of the point (ro , u) in the new coordinates 
with respect to the nontrivial (Ising) fixed point [i* which we are now approaching (see Fig. [3]). 
The singular part of the free energy density, / s , is after the transformation <f> — > ip denoted by / s , 

Furthermore, we introduce the quantity f s (fo,u,h) = f s (ro,U,h). Because the total free energy is conserved along 
the renormalization trajectory, the singular part of the free energy density changes as 



f ( '^0_ J£_ M _l-df ( I2_ ;2 _^_,4-d K\+d/2 

/s U 2 ' R 4 ' RJ \R 2 ' R 4 'R 

= ij-«/(4-d) / g (r R 2d /( 4 - d \u, hR 3d /V- d ^ 

= b -d R -M/(4-d)f s ^ R 2d/(4-d) b y t ^ b y^ hR M/(4-d) b y^ ^ (g) 



1 Naturally, this rescaling is not compulsory and the same results will be obtained without it, provided one keeps track of the 
dependence of the nontrivial fixed point on cR 2 , arising from the integration over propagators in the inner part of the Brillouin 
zone. 
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where we have used the notation t = [T — T C (R)]/T C (R) for fo- In Fig. pi t stands for the distance to the critical 
line connecting [1q and /i*. In the second equality we have substituted the crossover scale, I = R 4 ^ 4 - d \ Of course, 
this is only a qualitative measure for the location of the crossover, but the renormalization predictions for the scaling 
exponents are exact. The relation (^), which holds for 1 < d < 4, is the key to the scaling relations obtained 
on phenomenological grounds in Ref. B. We will first illustrate this by deriving the R dependence of the critical 
amplitudes of the magnetization density m and the magnetic susceptibility \. The magnetization density can be 
calculated by taking the first derivative of the free energy density with respect to the magnetic scaling field h, 

df S/ df s ( r u h 

= b V*-d R -d/{4-d)f(l) ^ R 2d/(4-d) b y t ^ b y %hR 3d/(4-d) b y^ ^ (1Q) 

where denotes the first derivative of / s with respect to its third argument. Here we have neglected any nonlinear 
relation between the magnetic scaling field h and the physical magnetic field. Furthermore, we have omitted a 
contribution from the analytic part of the free energy, because h only couples to the k = (uniform) mode. To 
extract the dependence of m on t and R from Eq. (fLOl) , we choose the rescaling factor b such that the first argument 
of the derivative of / s is equal to 1, i.e., b — i _1 /»t_R- 2d /bt(4-<i)] , and set the irrelevant variable u and the magnetic 
scaling field h equal to zero, 

m = 1*R<W->Q/Wfp (1,0,0) . (11) 

This result agrees with Ref. Eq. (34)]. In the same way we can calculate the magnetic susceptibility from / s by 
taking the second derivative with respect to h, 

d 2 f s , ,s d 2 f s ( r u h 



= b *y*-d R 2d/(4-d) f(2) (t R 2d/(i-d) b y t ~ b y^ hR 3d/(,4-d) b y^ (12) 

Choosing the arguments of the second derivative of / s as in Eq. ([[l]), we find 

X = f-V i? 2d(l- 7 )/(4-d)/ s (2) (lj Qj Q) t (13) 



in agreement with Ref. J3, Eq. (39)]. In Eqs. (O) and (13) we have used the well-known relations between the 



renormalization exponents and the critical exponents (see, e.g, Table 18 in Ref. Jl 1|) . 

The finite-size scaling behavior of thermodynamic functions can also be derived from this renormalization scenario 
by including a finite-size field 1/L as an additional argument of the free energy density in Eq. (^|). Under the first 
renormalization transformation this scaling field will scale as l/L and under the second renormalization transformation 
as Ib/L = 6i? 4 /^ 4_d ^i _1 . The finite-size scaling behavior is found by choosing b such that Ib/L = 1, i.e., b = 
LR-^ {i - d \ Substituting this into Eqs. © and (|l|), we obtain for m 



m = L V b -d R (3d-4y h )/(4-d) f(X) ^ R -2(2y t -d)/(4-d) R -4 w /(4-d) h[ y* R (3d-4y b )/(4-d) J ^ (u) 

and for \i 

x = tfy h -d R 2{3d-4y h )/{4-d) j(2) (+ L y t R -2{2y t -d)/{4-d) ^ uL Vi R~ 4yi/ {4 ~ d] , /iL^^(3d-4 ah )/(4-d)^j _ ^ 

These results agree with Ref. M, where the prefactors of the magnetization density and the magnetic susceptibility 
were predicted as, respectively, £-£/^(2/3-7)/M4-d)] and L1 /u R (4p -2^/^(4- d)\^ Furthermore, the first argument 
of the scaling functions was predicted as tL 1 /"R' c , with JC = -(2a)/[u(4 - d)] (Ref. g Eq. (25)]). This is indeed 
equivalent with our result K, = —2(2y t ~d) / (A—d). However, the predicted range dependence of the critical amplitudes 
(i.e., of the prefactors in Eqs. ( p"4| ) and (p"5|)) is only valid in the limit of infinite range. For smaller ranges, i?-dependent 
correction terms are present. These correction terms can be calculated as well. They do not come from the dependence 
of the scaling functions on the irrelevant fields, as corrections to scaling normally do: these corrections vanish in the 
thermodynamic limit. However, they come from higher-order contributions to the renormalization of the ip 4 coefficient 
which were previously neglected in the derivation of the crossover scale lo. Note that in the neighborhood of the 
Gaussian fixed point, the terms ip n with n < 2d/(d — 2) are relevant and that for d = 2 all higher-order terms are 
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equally relevant. However, the coefficients of these terms are, after the rescaling 4> — ► ip, proportional to R~ n , so the 
leading contribution comes from the term wR~ 6 ip e . Under a spatial rescaling with a factor I = e s , the renormalization 
equation for this term is, to leading order, 



d 4- = (6 - 2d)w' 
as 



(16) 



The solution of this equation, w'(s) — tue' 6 2d > s can be substituted in the renormalization equation for the - 4 
coefficient, 



f du' u' w' 



To first order in w, this yields, 



= JA-d)s J_ 



i? 4 



I? 4 



+ 2 - d R 2 V V 



i? 4 



2-dR? 



l 2 - d - 1) 



(17) 



(18) 



where it and w denote the values of u' and w' at Z = 1. 
crossover scale Zo = i? 4 /( 4 ~ d ) is multiplied by a factor (1 -j 
correction. However, the solution (|l^) is not valid for d = 
solution of Eq. (J17|) is then given by 



respectively. This implies that the previously obtained 
dR~ 2 ) and hence all critical amplitudes will exhibit this 
2, where Mi? _4 -0 4 and wR~ 6 ip 6 are equally relevant. The 



U> 5, 1 / W \ , 2 1 / W . 

R^ e R-*{ u+a & s ) = l ^r +a i^ ln 



(19) 



which yields a (leading) correction factor [1 + R~ 2 (di + &2 lni?)] in the crossover scale and the critical amplitudes. 

From a similar mechanism we can derive the R dependence of the so-called shift of the critical temperature ^, 
Eq. (15)]. A detailed treatment of the shift of T c can be found in, e.g., Refs. (l^||. It arises from the u-dependent 
contribution in the renormalization equation for the i/j 2 term, 



R? ds R 2+ R±- 
Thus, the first argument on the RHS of the first equation in (0) should be replaced by 



(20) 



'0 _ i2 

R 2 



fro_ _ ~uo_\ 
\R 2 °R i ) 



A 2 



1 

R 2 



To 



C R^) +C R 2 



(21) 



The term between round brackets is proportional to the reduced temperature and the last term is the leading shift. 
Substitution of the crossover scale Iq shows that this shift in the reduced temperature is proportional to R^ 2d /^ d \ 
which indeed vanishes in the mean- field (R — ► 00) limit. Remarkably, this disagrees with Ref. [Q, where a shift oc R~ d 
was predicted. Unfortunately, it is not possible to settle this issue at present by means of Monte Carlo data, because 
only results for d = 2 are available, which is a special case. Viz., for d = 2 we obtain instead of Eq. (|2|) the following 
solution of Eq. (EfJ) , 



'0 
R 2 



= I- 



\R 2 



c T 2lnl 



(22) 



Thus, we find, upon substitution of the crossover scale, that the shift in the reduced temperature has the form 
(p + oln R)/R 2 , where the constant p comes from a multiplicative factor introduced by the crossover criterion. In 
Ref. B, d — 2 was already suggested as a special case, with possibly logarithmic corrections. The renormalization 
argument shows that such a lni? term is indeed present. 

Now, let us return to Eq. (|]), where we omitted quartic (and higher) terms in kR. It follows from the renormalization 
scenario that terms proportional to k 2n transform as k 2n l 2 ~ 2n under the first renormalization transformation and hence 
are irrelevant for n > 1. The behavior of these terms under the second renormalization transformation is less simple, 
but again quartic and higher terms do not influence the leading terms, see, e.g., Ref. Jl3|, Section VII. 6]. 

Besides, it can be seen that the structure of th e in teraction term does not depend on the details of the spin-spin 
interaction. E.g., replacing the interaction term (Al) with K(r) = cR~ d exp [— (r/i?) 2 ] leads to precisely the same 
structure of the LGW Hamiltonian and hence to the same scaling relations involving R. This is in agreement with 
the universality hypothesis. 
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Furthermore, the renormalization description explains why the interaction range R must be large to observe the 
predicted powers of R in the critical amplitudes: only for systems with R large the renormalization trajectory starts in 
the neighborhood of the Gaussian fixed point and hence only these systems will accurately display the corresponding 
R dependence. Finally, in the finite-size scaling description, the system size must be sufficiently large in order to 
observe the crossover to Ising-like critical behavior: we require that the rescaling factor b is minimal of order unity, 
or L = e>(i? 4 /( 4 -«)- 



III. MONTE CARLO RESULTS AND COMPARISON WITH THE THEORETICAL PREDICTIONS 



A. Definition of the model 



We have carried out Monte Carlo simulations for two-dimensional Ising systems consisting of L x L lattice sites with 
periodic boundary conditions and an extended range of interaction. Each spin interacts equally with its z neighbors 



lying within a distance R m , as defined in Eqs. (|l|) and (Al) with R replaced by R m and d = 2. The Monte Carlo 
simulations were carried out using a special cluster algorithm for long-range interactions 0. Its application to the 
interactions defined above is described in Appendix |b|. Following Ref. Q we define the effective range of interaction 
R as 

r2 = E^-r,) 2 ^ (23) 
' |r, - r 3 1 2 with \n - Tj] < Rm . (24) 



z 71 



Table | lists the values of R 2 n for which we have carried out simulations, as well as the corresponding values of R 2 . 
The ratio between R 2 and R 2 ^ approaches 1/2, as can be simply found when the sums in Eq. (p3| ) are replaced by 
integrals. Note that the results for R 2 ^ = 18 and R 2 n = 32 cannot be compared to those presented by Mon and 
Binder, because in Ref. pj the interactions were for these two system sizes spatially distributed within a square, as 
can be seen from the number of interacting neighbors and the corresponding effective ranges of interaction. 



B. Determination of the critical temperature 



The critical temperatures T c of these systems have been determined using the well-known universal amplitude ratio 
Ql = ( m2 )i/( TO4 ) L- Both in the Ising and in the mean-field limit the critical-point value of this quantity is accurately 
known. In the mean-field limiiLQ = Qmf = 87r 2 /[r(±)] 4 « 0.456947, see Refs. |l2) and [0, Appendix A]. In the Ising 
limit, Q = Qi = 0.856216 (1) Q. 

As was noted in Section || and also in Ref. [Q, rather large system sizes (0(R 2 )) are required to determine T c , 
since Q must approach Q\. For R 2 n < 10 we have included linear system sizes up to L = 500 and for larger ranges we 
have used system sizes up to L — 700 or even L = 800 (i?^ = 100, 140). For each run we have generated 10 6 Wolff 
clusters after equilibration of the system. The various thermodynamic quantities were sampled after every 10th Wolff 
cluster. In Fig. g, Ql(K c ) for i?^ = 140 is plotted as function of the system size. One clearly observes the crossover 
from Qmf (for L < R 2 m ) to Qi (for L > R^). 

In the Ising limit, the finite-size expansion of Ql reads (see, e.g., Ref. |ll|] ) 

Q L {K) = Q + ax(K - K c )L m + a 2 (K - K C ) 2 L 2 ^ +■■■ + b x L yi + b 2 L y * +■■■ . (25) 

K denotes the spin-spin coupling, K c the critical coupling, and the dj and bi are non-universal (range-dependent) 
coefficients. The term proportional to L V2 , with j/2 = d — 2yh, comes from the field dependence of the analytic 
part of the free energy. In a 4 theory this term is absent, as was stated in section ||, but in a discrete model it 
should be allowed for. The exponents yt, yh, and ?/;, which have already been introduced in the previous section, are, 
respectively, the temperature, magnetic, and leading irrelevant exponent for the two-dimensional Ising model; yt = 1 , 
2/h = and j/i = —2. Table || displays the results of a least-squares fit according to Eq. (|5|), where y t , yi, and 
?/2 were kept fixed at their theoretical values. For comparison we have included the estimates for K c from Ref. [Q. 
Except for i? 2 ^ = 10 there is good agreement between the respective estimates. The discrepancy for i?^ = 10 may 
be explained by the limited range of system sizes in Ref. Q. Furthermore, for R^ = 2, which corresponds to the 
Ising model with nearest and next-nearest neighbor interactions, an accurate transfer-matrix estimate of the critical 
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coupling exists, K c = 0.19019269(5) [Gl|. The Monte Carlo result agrees with this value. The results for Q are in 
good agreement with the expected value Qi, which confirms not only that universality is satisfied, but also that the 
maximum system sizes in our simulations are sufficiently large, so that crossover to Ising-like critical behavior indeed 
has taken place, as it should for an accurate determination of K c . In fact, the error margins on K c can be reduced 
significantly by fixing Q at its Ising value in Eq. (pq ) (see Table [n]). Figure ^ illustrates the shift in l/(zK c ) oc T c 
as function of R~ 2 . Even close the mean-field limit (R~ 2 — > 0), the deviation of l/(zK c ) from 1 seems not truly 
linear. Therefore we have tried to identify the logarithmic term, which was suggested in Rcf. [Q and derived from the 
renormalization scenario, by writing the following expression for the critical coupling, 

p + q\nR 

zK c = l + ■ (26) 

In Fig. H we have plotted A = (zK c — \)R 2 versus lni?. Indeed, for large values of R the points lie approximately on 
a straight line, confirming the presence of the logarithmic correction. 

Another In R correction was suspected in Ref. B in the temperature-dependent argument of the finite-size scaling 
functions. This argument is proportional to R-^vt-d)/(i-d) = R -2 a /[»(4-d)] ^ see Eqg _ Q and For = 2, 

a = implies a logarithmic divergence of the specific heat and hence one might expect a similar logarithmic factor 
here. On the other hand, we have not found a mechanism in the renormalization scenario which could explain such 
a factor. Therefore, we have numerically examined the range dependence of the coefficient a\ in Eq. (|2f|). Since 
(K — K c ) is proportional to R~ 2 , we must first divide a\ by R 2 . Figure |5| displays this quantity as function of the 
range. For small ranges, there is a strong dependence on R, but the coefficients seem to approach a constant value in 
the large-range limit. This suggests that a logarithmic correction factor is absent. 



C. Range dependence of the magnetization density 

We have sampled the absolute magnetization density, (|m|), for which the range dependence is given by Eq. (fh|). 
This quantity has been fitted to the following finite-size expansion, 

m L {K, R) = L y *~ d {d (R) + dx(R)[K - K c (R)]L Vt + d 2 (R)[K - K c {R)] 2 L 2yt + ■■■ + ei (R)L v ^ + • ■ ■} , (27) 

where we now have explicitly indicated the range dependence of the parameters. The critical couplings found from 



this quantity agree well with those obtained from the amplitude ratio Q and the exponent y^, listed in Table HI, is in 



good agreement with the exact value 15/8. Furthermore, we have made a least-squares fit with K c fixed at the most 



accurate values obtained from Q. The corresponding estimates for j/h are also shown in Table III. They lie even closer 
to 15/8, which again corroborates that all systems belong to the Ising universality class. From the critical amplitudes 
do(R) we can derive the leading R dependence of the magnetization. To increase the accuracy, the values in Table III 
were determined with fixed at its theoretical value. As can be seen from the log-log plot in Fig. ^|, the approach 
to the asymptotic scaling behavior is very slow. Therefore we have determined the scaling exponent in two different 
ways. A straight line through the points for the three largest ranges yielded do(R) oc ^-°- 738 ( 13 ) ; i n agreement with 
the predicted exponent —3/4 (Eq. (pT[)). Inclusion of the correction factor [1 + R~ 2 (di + a 2 hii?)], as predicted from 
Eq. (|l9|), allowed us to include all data points in the fit and yielded dg{R) oc /j>-°- 756 (5) ^ a j so m g 00C j agreement with 
the predicted exponent. 



D. Range dependence of the susceptibility 

The magnetic susceptibility can be calculated from the average square magnetization, 

X = L d (m 2 ) . (28) 

We thus expect the following finite-size scaling behavior, 

Xl(K,R) = s Q + L 2 y*- d {p (R)+ Pl (R)[K - if c (i?)]i yt + P 2(R)[K - K C (R)] 2 L^ + ■■■ + qi(R)L* + • • •} . (29) 

The term sq comes from the analytic part of the free energy. Because it tends to interfere with the term proportional 
to qi(R), we have ignored it in the further analysis. Again, the critical couplings obtained from a least-squares fit lie 
close to those in Table || and the estimates for y^ agree with the Ising value (see Table |Tv| ). By repeating the fits with 
K c fixed at the most accurately known values, the values for y^ lie even closer to 15/8 (third column of Table |rv|). 
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From the parameter po(R), plotted in Fig. fj], we can extract the leading range dependence of the susceptibility. A 
straight line through the amplitudes for the three largest ranges gave Po(R) oc i? -146 ^ 3 ^. For a curve (including the 
first correction term) through the amplitudes it was necessary to include the data for all ranges R 2 > 7/3 in the fit, 
in order to determine the coefficient of the In R factor. This yielded po(R) oc Rr 1A7ty2 K Both exponents are in good 
agreement with the predicted value 2 (3d — %h)/(4 — d) = —3/2. 



E. Spin— spin correlation function 

The finite-size scaling behavior of the spin-spin correlation function g(r) closely resembles that of the magnetic 
susceptibility \, as may be expected from the fact that x is the spatial integral of g. We also expect the range 
dependence of the two quantities to be the same. We have sampled the correlation function over half the system size 
and analyzed it using the expansion 

g L (K, R) = L 2 y*- 2d {v (R) + v x {R)[K - K C {R)]L^ + v 2 {R)[K - K C (R)] 2 L^ +■■■ + wx(R)U* + • • •} . (30) 

The constant term in ( p9f ) is not present here (see, e.g, Ref. (Tl]]). Table shows the results for yh, both with K c 
free and fixed. In the latter case, j/h is in accurate agreement with its theoretical value, just as for the magnetization 
density and the magnetic susceptibility. Figure || shows a log-log plot of the critical amplitude v (R) as function of 
the range. A fit of a straight line through the points with R 2 > 35 (i.e., R^ > 72) yielded vo{R) oc R~ 1A& ( 3 \ whereas 
a curve through all points with R 2 > 7/3 gave Vq(R) oc R~ 1A9 ( 2 \ Both estimates are again in good agreement with 
the predicted exponent —3/2. 



IV. CONCLUSION 



In this paper, we have derived the dependence of scaling functions on the range of interactions from renormalization- 
group arguments. The results agree with the predictions of Mon and Binder and yield in addition the corrections 
to the leading scaling behavior, including the previously conjectured logarithmic factor in the shift of the critical 
temperature of two-dimensional systems. 

We have also carried out accurate Monte Carlo simulations for systems in which the range of the interactions was 
large enough to verify the theoretical predictions. It was confirmed with high precision that all examined systems 
belong to the 2D Ising universality class. Besides the range dependence of critical amplitudes, we also observed the 
predicted range dependence of the corrections to scaling. 



APPENDIX A: FOURIER TRANSFORM OF A SPHERICALLY SHAPED INTERACTION PROFILE 

We define the following isotropic spin-spin interaction (the subscript d denotes the dimensionality), 

KJr) = { cR ~ d if I 1 "' ^ R (Al) 
Kd(r) ~ \ if |r| > R (AL) 

We have normalized the interaction strength, to make the integrated interaction (energy) independent of the range. 
In this appendix, we calculate the Fourier transform of this interaction for a general number of dimensions. For d = I 
the calculation is trivial: 

r r +R 2r 
Ki(jfc) = - J dx e ikx = — sm(kR) . (A2) 



For d = 2 and d — 3 one obtains Bessel functions. Using the equality J\/2{x) = y/2/(Trx) sin(a;), the results for 
d = 1,2,3 can be summarized as 



Kd{k) = c \Jr) J ^ kW >' < A3 > 

where J„ is a Bessel function of the first kind of order v. This suggests that this equality is valid for general d, which 
can indeed be shown by induction. If we assign the x coordinate to the (d + l)th spatial dimension and use the 
notation k 2 = Xw=i ^li we m ay write 
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n+R / n \ d/2 

K d+1 (k) = J dxcos(k x x) (R 2 - x 2 fl\j dl2 (kdVR 2 -* 2 ) 

2c [2-K\ d/2 f R / \ p( d + 2 )/ 2 



2jT \(d+l)/2 



= cf — j J(d+i)/2(kR) , (A4) 

where we have used a Hankel transform of general order, see, e.g., Ref. [pl| p. 40, Eq. (48)]. 



APPENDIX B: MONTE CARLO ALGORITHM FOR SPIN SYSTEMS WITH MEDIUM-RANGE 

INTERACTIONS 

The cluster algorithm we have used for the present Monte Carlo simulations has essentially been described in 
Ref. @. However, that description is rather concise; here, we present a somewhat more elaborate discussion of the 
mathematical backgrounds, and we outline how the algorithm is applied to medium-range models of the type studied 
in this paper. 

The description is given in terms of the Wolff cluster algorithm |l7| ], but the principle applies only to the cluster 
formation process. Thus, it is also applicable in the Swendsen-Wang case jl8f| . For simplicity, we describe the way a 
cluster of spins is built in the case that there are only ferromagnetic interactions. 

For each spin in the cluster, we have to run a task described below. During this task, new spins may be included in 
the cluster. For this reason, it is convenient to use a "stack" memory containing the addresses of the spins for which 
the task remains to be done. 

The task for a spin Si read from the stack is the following. A loop is executed over all neighbors sj interacting with 
Si. In each step of this loop, the bond connecting sites i and j is "activated" with a probability 

p(s l ,s j ) = S SlSj p , (Bl) 

where p = 1 — exp(— 2K), in which K is the coupling between Si and Sj. The simulation process would conventionally 
include a test whether Si and Sj are parallel, and if so, the production of a uniformly distributed pseudo-random 
number r. If r < p, the bond is activated, Sj is added to the cluster, and its address j is stored in the stack memory. 
Since this loop runs over all neighbors interacting with Si, i.e., over all sites within a distance R m , the process 
becomes very time consuming when the range R m of the interactions becomes large, just as in the case of Metropolis 
simulations. 

However, the cluster formation process can be formulated in a more efficient way. Part of the work involved in the 
activation of the bonds between Si and its neighbors Sj can be done in a way that does not depend on the signs of the 
spins. Thus, as a first step, the bonds connected to Si are "provisionally activated" with a probability p, independent 
of their relative sign. Typically, only a small number of bonds will be provisionally activated for each entry i in the 
stack memory (i.e., each spin in the cluster). Then, in the second step, the provisionally active bonds, say between 
sites i and j, are actually activated if = Sj , i.e., with probability S SiS . . During the second step, the bonds that were 
left inactive in the first step can be ignored. 

Since the first step does not depend on the signs of the spins, and the probability p is typically quite small, we 
introduce (following Ref. JtJ) a cumulative bond probability. This quantity determines which bond is the next bond to 
be provisionally activated. The probability that, during the first step, k — 1 bonds are left inactive and that the kth 
bond is provisionally activated is equal to 

P(k) = (l-p) k - 1 p. (B2) 
The cumulative bond probability C(k) is defined as 

k 

C(k) = J2 P(n) ■ (B3) 

n=l 

The interval fc — 1 to the next bond to be provisionally activated is obtained by drawing a random number r (g [0, 1)). 
If this random number lies between C(fc — I) and C(fc), k — 1 bonds are skipped and the fcth bond is provisionally 



activated. It is readily seen that this procedure leads to precisely the required probabilities given in Eq. (B2). But 
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the number of operations per spin in the cluster is only of order pR^l near criticality, this quantity is approximately 
equal to 1. Thus, the work involved in the decision concerning the actual bond activation is also of order unity. 

We check independently whether the resulting probability of activating the first bond at a distance k is equal to 
that in the conventional approach. Consider a cluster spin, say So, with a chain of neighbors denoted si . . . Sk, of 
which to spins are antiparallel to sq and I = k — m spins are parallel to sq, among which Sk- In the conventional Wolff 
cluster algorithm, the probability that Sk is the first spin to be added to the cluster is, provided that sq — s&, given 
by 

l m (l-p) ,_1 p=(l-p) , - 1 p. (B4) 



On the other hand, if we use the cumulative bond probability (B3), this probability is calculated as follows. Either the 
fcth spin is selected directly (if the first random number lies between C(k — 1) and C(k)) or one of the to antiparallel 
spins is selected, say s a , which is of course not added to the cluster. In the latter case, another random number is 
drawn and a new spin is selected. Again, this may be the fcth spin, or one of the remaining antiparallel spins between 
s a and Sfc. Now, let us show that the sum of these probabilities of adding Sk as the first spin to the cluster is equal 



to (B4). Denote the number of selected, "intermediate", antiparallel spins by i. There are (" l ) possibilities of selecting 
i intermediate spins. The probability of selecting Sk after each of these sequences of spins is p l (l — p) m ~ l (l — p) l ~ 1 p. 
The total probability is the sum over all numbers of intermediate spins 

m 

£ (T) - f) m_< (l - P) ,_1 P = (1 -P)^P > (B5) 

i=0 



which is indeed equal to (B4) 



As shown in Ref. 0, Cj(k) = 1 — exp(— Yln=j+i ^0 = ^ — cxp(— [A; — j]K). By inverting this relation, the bond 
distance k can be calculated from Cj(k), i.e., from the random number. This approach is highly efficient; compared to 
conventional (Metropolis) algorithms the gain is typically a factor 0(Rf n L 2 ). Finally we remark that efficient variants 
of this technique can be applied to long-range 0(n) models with n > 1. Again, bonds are provisionally activated 
with a probability p = 1 — exp(— 2K); actual activation of a bond between spins st and sj is done afterwards with 
a probability [1 — exp(— 2Ksi iX 8j tZ )]/p if (s*$ • Sj) > (where z defines the spin-flip direction of the cluster-formation 
step) and otherwise with probability 0. 
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FIG. 1. Qualitative picture of the renormalization trajectory describing the crossover from the Gaussian fixed point 
fio — (0, 0) to the Ising fixed point (/,* — (r\ h u*). 




1000 



FIG. 2. The critical-point amplitude ratio Ql(K c ) for = 140 as function of the system size. For large L, Ql{K c ) 
approaches the Ising limit Qi ~ 0.856216. For small L, Ql(K c ) approaches the mean-field limit Qmf ~ 0.456947, although 
there are significant finite-size corrections for the smallest values of L. 
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0.02 0.04 0.06 0.08 0.1 

R- 2 

FIG. 3. Plot of l/(zK c ) versus R~ 2 . The dashed line denotes the extrapolation to the mean-field limit. The inset shows 
l/(zK c ) over the full range of R~ 2 between the Ising and the mean-field limit. 
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FIG. 4. A = (zK c — l)R 2 versus In 7?. For large R the graph strongly suggests the presence of a logarithmic correction in 
the shift of the critical temperature. The error bars do not exceed the symbol size. 
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FIG. 5. Range dependence of the amplitude of the temperature-dependent argument of the finite-size scaling function of the 
universal amplitude ratio Q. 
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FIG. 6. The critical amplitude do(R) of the magnetization density versus R 2 
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FIG. 7. The critical amplitude po(R) of the magnetic suscetibility versus R 2 
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FIG. 8. The critical amplitude vo{R) of the spin-spin correlation function versus R 2 
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JTL m 


z 
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3 
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20 


17 

5 
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25 
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224 
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4003 
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TABLE I. The range of interaction R m , the corresponding number of neighbors z and effective range of interaction R. 



r>2 


Q 


K c 


K c 


Kc p 


2 


0.8556 (5) 


0.1901908 (19) 


0.1901931 (11) 


0.190 


4 


0.8557 (9) 


0.1140216 (18) 


0.1140225 (7) 


0.11402 


6 


0.8553 (7) 


0.0631917 (8) 


0.0631926 (4) 




8 


0.8553 (13) 


0.0510460 (10) 


0.0510467 (4) 


0.05106 


10 


0.8563 (9) 


0.0324136 (5) 


0.03241352 (18) 


0.032463 


18 


0.8555 (14) 


0.0185335 (3) 


0.01853367 (9) 




32 


0.853 (3) 


0.01075152 (25) 


0.01075182 (7) 




50 


0.856 (6) 


0.00657274 (26) 


0.00657276 (5) 




72 


0.854 (4) 


0.00464056 (16) 


0.00464064 (4) 




100 


0.850 (8) 


0.00325903 (15) 


0.00325905 (5) 




140 


0.862 (17) 


0.00234637 (19) 


0.00234631 (2) 





TABLE II. The amplitude ratio Q and critical coupling K c for the various ranges of interaction studied in this paper. The 
numbers in parentheses denote the errors in the last decimal places. The third column shows the estimates for K c obtained 
with Q fixed at its Ising value. For comparison, we also list the estimates for K c given in Ref. W- 



r>2 


Vh 


Vh 


d (R) 


2 


1.8745 (7) 


1.8749 (3) 


0.9533 (4) 


4 


1.8763 (15) 


1.8756 (4) 


0.8706 (5) 


6 


1.873 (3) 


1.8767 (13) 


0.7937 (10) 


8 


1.873 (3) 


1.8754 (8) 


0.7523 (7) 


10 


1.874 (2) 


1.8748 (7) 


0.6783 (6) 


18 


1.871 (3) 


1.8740 (12) 


0.5816 (6) 


32 


1.875 (6) 


1.8744 (9) 


0.4929 (11) 


50 


1.873 (7) 


1.876 (2) 


0.4181 (18) 


72 


1.865 (5) 


1.8752 (16) 


0.3742 (8) 


100 


1.867 (9) 


1.877 (2) 


0.3296 (8) 


140 


1.895 (13) 


1.879 (3) 


0.2938 (13) 



TABLE III. The exponent yh and the critical amplitude of the magnetization do(R) for the various ranges of interaction. 
The third column shows the estimates for yh obtained with K c fixed at the most accurate values shown in Table n\. 



15 



d2 


Vh 


Vh 


Po(R) 


2 


1.8754 (9) 


1.8748 (2) 


0.9743 (9) 


4 


1.8753 (12) 


1.8752 (3) 


0.8136 (7) 


6 


1.8740 (18) 


1.8761 (10) 


0.6762 (14) 


8 


1.873 (2) 


1.8750 (6) 


0.6076 (9) 


10 


1.874 (3) 


1.8741 (6) 


0.4943 (7) 


18 


1.874 (4) 


1.8740 (11) 


0.3620 (9) 


32 


1.868 (4) 


1.873 (2) 


0.2622 (9) 


50 


1.862 (6) 


1.874 (3) 


0.1914 (7) 


72 


1.863 (17) 


1.870 (4) 


0.1534 (8) 


100 


1.870 (6) 


1.874 (4) 


0.1180 (8) 


140 


1.86 (3) 


1.870 (5) 


0.0954 (9) 



TABLE IV. The exponent j/h and the critical amplitude of the magnetic susceptibility po(R) for the various ranges of 
interaction. The third column shows the estimates for yh obtained with K c fixed at the most accurate values shown in Table |ll 



d2 


Vh 


yh 


v (R) 


2 


1.8759 (8) 


1.8754 (3) 


0.9405 (9) 


4 


1.8744 (12) 


1.8750 (3) 


0.7860 (6) 


6 


1.8748 (19) 


1.8765 (11) 


0.6528 (15) 


8 


1.8746 (17) 


1.8754 (6) 


0.5862 (9) 


10 


1.875 (3) 


1.8741 (7) 


0.4770 (9) 


18 


1.874 (4) 


1.8745 (10) 


0.3489 (9) 


32 


1.873 (4) 


1.8747 (18) 


0.2541 (6) 


50 


1.864 (7) 


1.876 (3) 


0.1827 (9) 


72 


1.860 (8) 


1.871 (4) 


0.1473 (9) 


100 


1.872 (9) 


1.876 (4) 


0.1128 (9) 


140 


1.86 (3) 


1.871 (3) 


0.0915 (9) 



TABLE V. The exponent j/h and the critical amplitude of the spin-spin correlation function vo(R) for the various ranges of 
interaction. The third column shows the estimates for y^ obtained with K c fixed at the most accurate values shown in Table ffR 
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